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ABSTRACT 

The discrepancy between abundances computed using optical recombination lines (ORLs) 
and collisionally excited lines (CELs) is a major unresolved problem in nebular astrophysics. We 
show here that the largest abundance discrepancies are reached in planetary nebulae with close 
binary central stars. This is illustrated by deep spectroscopy of three nebulae with a post common- 
envelope (CE) binary star. Abell 46 and Ou5 have abundance discrepancy factors larger 

than 50, and as high as 300 in the inner regions of Abell 46. Abell 63 has a smaller discrepancy 
factor around 10, but still above the typical values in ionized nebulae. Our spectroscopic analysis 
supports previous conclusions that, in addition to “standard” hot (Te'^lO'^ K) gas, a colder 
(Te^lO^ K) ionized component that is highly enriched in heavy elements also exists. These 
nebulae have low ionized masses, between 10“^ and 10“^ Mq depending on the adopted electron 
densities and temperatures. Since the much more massive red-giant envelope is expected to be 
entirely ejected in the CE phase, the currently observed nebulae would be produced much later, 
in post-CE mass loss episodes when the envelope has already dispersed. These observations add 
constraints to the abundance discrepancy problem. Possible explanations are revised. Some are 
naturally linked to binarity, such as for instance high-metallicity nova ejecta, but it is difficult 
at this stage to depict an evolutionary scenario consistent with all the observed properties. The 
hypothesis that these nebulae are the result of tidal destruction, accretion and ejection of Jupiter- 
like planets is also introduced. 

Subject headings: planetary nebulae: individual (A 46, A 63, Ou5) - ISM: abundances - binaries: close 
- novae, cataclysmic variables - planet-star interactions 


1. Introduction 


This work deals with two apparently unrelated, 
main topics in the study of planetary nebulae 
(PNe). 

The first one is the role of binary evolution, 
which is the favoured explanation of the diverse 
morphologies display ed by PNe (e.g. ISoker 1997 : 
Balick fc Frank 2002fl . but whi ch may even be the 


cause of their me re existence (|Moe fc De Marco 
20011: ISokeilliinbll . 


The second issue is the so-called abundance dis¬ 
crepancy problem. It is well known (see e.g. Os- 
terbrock & Ferland 2006) that in photoionized 
nebulae - both PNe and H II regions - optical 


recombination lines (ORLs) provide abundance 
values that are systematically larger than those 
obtained using collisionally excited lines (CELs). 
This a long-standing problem in nebular astro¬ 
physics, and has obvious implications on the mea¬ 
surement of the chemical content of the Universe, 
often done using CELs from emission regions and 
the ISM. The abundance discrepancy factor {adf) 
between ORL s and CELs is usually between 1.5 
and 3 (see e.g. Garcia-Roias fc Esteban|[2007t iLiul 


2012 : Esteban et al.ll2014l) . but in PNe it has a sig¬ 

nificant tail extending to much larger v alues. The 
nebu la with the largest known adf (~70,lLiu_et_aL 
20061) is Hf 2-2, whic h has a close binary central 
star ( Lutz et ffi]ll998ll . 
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In our recent study of the new Galactic PN 
IPHASXJ211420.0+434136 (Ou5~), which also 


has a binary central star ( Corradi et al. 2014h . 


we noticed the unusual strength of recombina¬ 
tion lines such as C Ii4267 A and the 4650 A 
O IH-N III blend. This points to large adf val¬ 
ues also in this object, as well as i n other binary 
PNe such as Abell 46 (Bondi lOSoll and Abell 63 


( Bond. Tiller fc Mannerv 1978ll in which stron g 


C II emission was detected ( Pollacco fc Belll[l9^ ). 
The short orbital periods, between 8 and 11 hours, 
indicate that the central stars of these PNe have 
gone through a common-envelope (CE) phase in 
their previous evolution. Therefore, we decided to 
obtain deep spectroscopy of these three nebulae 
to investigate whether the abundance discrepancy 
problem in PNe is in some way related to the bi¬ 
nary nature of their central stars. That this is 
indeed true is demonstrated by the results pre¬ 
sented in the following. 


2. Observations 


Ha-|-N II images of A 46 and A 63 were retrieved 
from the archives of the Isaac Newton Group of 
Telescopes and ESO. The image of A 46 is an 
1-hr exposure from the 2.5m Isaac Newton Tele¬ 
scope on La Palma obtained in 2008 with the Wide 
Field Camera (pixel size 0".33). Seeing was 1".0 
FWHM. The image of A 63 was obtained in 1995 
at the 3.5m ESO/NTT with EMMI (exposure time 
30 min, pixel size 0".27, seei ng 0".9). The image o f 
Ou5 is the one described in Corradi et al.l ( 2014ll . 
The inner regions of the nebulae are displayed in 
Figured] 


Spectra of A 46, A 63 and Ou5 were obtained 
on August 16-18 2014 with the 4.2m WHT tele¬ 
scope and the double-arm ISIS spectrograph. The 
long slit of ISIS was opened to 1" width and po¬ 
sitioned through the brightest parts of the nebu¬ 
lae as indicated in Figure [T] and in Table dl In 
the blue arm of ISIS, grating R1200B was used, 
providing a dispersion of 0.22 A pix“^ and a res¬ 
olution of 0.8 A. In the red arm, grating R316R 
gave a dispersion of 0.92 A pix“^, and a resolution 
of 3.4 A. Two different grating tilts were adopted 
in different nights in order to have a global wave¬ 
length coverage from 3610 to 5050 A in the blue, 
and from 5400 to 9150 A in the red, with little 
vignetting at the wavelength ends. The spatial 



Fig. 1.— Ha-I-N II images of the nebulae. The 
field of view is 90"x90" in all frames. The loca¬ 
tion of the long slit is shown. In A 46, the green 
cross indicates where emission of recombination 
lines peaks along the slit. The position of the LIS 
identified in the nebula is also indicated. 
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scale was 0"2 pix“^ in the blue, and 0"22 pix“^ 
in the red. Depending on the wavelength setting, 
total exposure times were of 2 hours for A 46, 2 
to 3.3 hours for A 63, and 2 to 4.7 hours for Ou5. 
The spectrophotometric standar ds BD+28 4211, 
BD+33 2642, and Feige 110 from lOke ( 199di . and 
HR 718 from Hamuv et al.l (1994li were observed 


for flux calibration. Seeing varied significantly 
during observations, and only data obtained with 
a seeing value <1'.'6 FWHM were retained. Images 
and spectra were reduced using packages ccdred 
and twodspec in IRAF0. 


3. Data analysis 


To increase the signal from these low surface 
brightness nebulae, 1-D spectra were extracted 
by integrating all the nebular emission included 
in the long slit. Note that in A 46, at the po¬ 
sition indicated by the arrow in Figure [TJ the 
slit intersects an elongated knot. In the spectra, 
this feature clearly stands out for its strong [N ii] 
6548,6583, [S il] 6716,6731, and more moderate 
[O il] 3726,3729 emission relative to the surround¬ 
ing nebula, from which it seems to be blue-shifted 
by ~10 km s“^. It could be a low-ionization 
small-scale structure (LIS) as often foun d in PNe 
( Goncalves. Corradi. &: Mamoasol 2001 1. Other 
similar blobs are seen in the image, but it should 
also be noted that the field is also rich in small 
background galaxies. This feature does not affect 
the following analysis, and is not considered any 
further in this work. 


Figure [5] shows an illustrative portion of the 
spectra of the nebulae, in a wavelength range 
where some of the most relevant ORLs and CELs 
are located. Emission-line fluxes were measured 
by means of multi-Gaussian fit using splot. Both in 
the blue and red, fluxes measured in the two spec¬ 
trograph settings were rescaled using the emis¬ 
sion lines in the overlapping spectral range. The 
flux differences in these overlapping lines (after 
scaling) were used to estimate the errors in the 
flux measurement. Identification of the observed 
lines, their fluxes and errors are presented in Ta¬ 
ble [2l All line fluxes are normalised to H/3=I00, 


where the observed H/3 flux integrated along the 
slit is 3.53xI0-l^ I.SOxIQ-^'^, and 1.37x10-^'^ 
erg cm“^ s“^ for A 46, A 63, and Ou5, respec¬ 
tively. 


The nebular reddening was computed from the 
hydrogen line ratios, specifically the Balmer Ha, 
H/3, Hy, and Hd lines and the Paschen 10-3, II- 
3 and 12-3 trans itions. The reddening law of 
Fitzoatri^ ( 2004 1 was adopted. The theoreti¬ 


cal hydrogen line ratios have a small but non- 
negligible dependence on the electron density and 
temperature rie and Tg. As discussed below, 
in addition to a standard nebular component at 
Te~I0'^ K, in these nebulae there is evidence for 
a much colder component at ssIO^ K, as for in¬ 
stance indicated by the Balmer jump in A 46 
(Sect. S]). We therefore computed the reddening 
by adopting either Te=I2500 K, characteristic of 
the GELs, or Te=I000 K, as suggested by the 
ORLs indicators. The mean values of the loga¬ 
rithmic extinction at H/3 of each nebula (c(H/3)) 
are listed in Table [3] for both temperature assump¬ 
tions. Those derived using the CELs temperature 
are 0.17 dex larger, and i n good agreement with 
the previous estim a tes by iPollacco fc Belli (lI997^ 
and ICorradi et al.l (2014), as expected consider¬ 
ing that they were determined under the same as¬ 
sumption. In the case of A 46, the y are also con¬ 
sisten t with the reddening maps of ISchlegel et al 


( 1998l l. which indicate a total foreground dust col¬ 
umn corresponding to c(H/3)=0.17. In the follow¬ 
ing, we conservatively use the c(H/3) values de¬ 
termined with the CELs temperature to deredden 
the observed line fluxes, as the ORLs temperature 
are more uncertain, especially for A 63 and Ou5. 
However, we also repeated the whole analysis with 
the lower c(H/3) values resulting from the ORLs 
temperature. Differences in the results are small 
and do not affect our conclusions. Errors on the 
dereddened fluxes includes the contribution of the 
uncertainty in the c(H/3) value. 


4. Physical and chemical properties of the 
nebulae 

4.1. Physical conditions 


^IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni¬ 
versities for Research in Astronomy (AURA) under coop¬ 
erative agreement with the National Science Foundation. 


The physical conditions for our nebulae were 
computed using several emission-line ratios. The 
electron temperatures, Tg, and electron densities, 
rie, are presented in Table [H The computa- 
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Fig. 2.— Illustrative portion of the spectra of the nebulae integrated along the slit. 


tions of physical condit ions were done with PyNeb 
(ILuridiana et al.ll2014l) . a python based package 
that allows derivation of the physical conditions 
and ionic and elemental abundances from emis¬ 
sion lines. The methodology followed for the 
derivation of r ip, and Tp has been described in 


Garcia-Roias et all ( 2012ll . We used the state-of- 
the-art atomic data listed in Table 21 Errors in 
these diagnostics were computed via Monte Carlo 
simulations. 


We computed ne([0 il]) and ne([S il]) for the 
three nebulae, and with lower precision ne([Ar iv]) 
for A 46 an d Ou5 . Similarly to what found 
by iLiu et aP ( 20061) for Hf 2-2, we find that 


ne([0 il]) is significantly higher than ne([S il]) and 
ne([Ar iv]). For the following analysis, we adopted 
the weighted average of ne([0 ii]) and ne([S ii]) 
for A 46 and A 63, and of ne([0 il]), ne([S ii]) and 
ne([Ar iv]) for Ou5. Note that the main results 
are almost independent on the assumed Tip. 
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Fig. 3.— Spectrum of A 46 at the Balmer jump. 


The intensity of the auroral [N il] A5755 line and 
the trans-auroral [O il] AA7320+30 lines were cor¬ 
rected for the contributi on of recombinat ion using 
equations (1) and (2) in iLiu et all (l20n(lll . A pre¬ 
liminary calculation of the N^+/H+ and 
abundance ratios by means of N ii multiplet 3 
and O II multiplet 1 was used to compute the re¬ 
combination contribution to [N ii]5755 and [O ii] 
7320-1-30, respectively. This resulted (Table |6]) 
in relatively small corrections for the [N ii] line 
flux, but very large corrections for the [O il] line 
fluxes. This means that Tg from the correspond¬ 
ing temperature diagnostics would be lower than 
computed neglecting recombination. In particu¬ 
lar, the fact that recombination dominates the 
trans-auroral [O ii] emission (100% in A 46 and 
Ou5, and 56% in A 63), no reliable Tg estima¬ 
tion can be derived using these lines. Tg([N ii]) is 
also very uncertain owing to the large errors in the 
measurement of the faint [N ii] 5755 and multiplet 
3 N II lines. 


Given the relatively high ionization degree of 
our three objects, we also computed the contribu¬ 
tion of recombination to t he auroral [ O iii| A4363 
line using equation (3) in Liu et ahl ( 2000h . To 
estimate the 0^+/H+ ratio, we assumed 


03+/H+ = (He/He+)2/3 x (0+/H+ -h 0^+/H+). 


to be between ~14% and ^43% depending on the 
nebula for the slit-integrated spectra (Table El but 
see also discussion in Section IT^ . The values ob¬ 
tained for Tg with and without recombination cor¬ 
rection are summarized in Table El Owing to the 
large uncertainties in Te([N il]) and Tg([0 il]), and 
considering the overall high excitation of the nebu¬ 
lae, we finally adopted Te{[0 ill]) as representative 
of the whole nebula, (see Table El and discussion 
below). 

To test whether in these nebulae part of the 
ORLs emission comes f rom a cold , high- metallicity 

e also 
ORLs. Un¬ 


gas component (e.g. iLiu et all l2006l l 
checked some Tg diagnostics from 
fortunately, the most widely used diagnostics, 
that is the ratio between the Balmer disconti¬ 
nuity at 3646 A and the Balmer decrement ratio, 
could only be measured in A 46 (Figure EJ- In 
this nebula, the Balmer continuum temperature 
in K) was dete rmined following the equation by 


Liu et al.l (|2001ll : 


/ Rpic 

r(Bac) = 368 x(l-h0.259?/+-f 3.40 V+) ■— 

V Hll 


- 3/2 


where ?/+ and are the He'''/H+ and He^'*'/H+ 
ratios respectively, and Rac/Hll is the ratio of 
the discontinuity of the Balmer jump in erg cm“^ 
s“^ to the Hi Hll line flux. For A 46, a very 
low Tg(Bac) = 1150 ±550 K, is obtained. We also 
considered other Tg line diagnostics from recombi¬ 
nation lines, such as Tg derived from He i line ra¬ 
tios (Tg(He I 5876/447 1, 6678/44711. rZhang et al 


20051 : iLiu et ah 2006ll and Tg derived from the 
ratio of O ii ORLs to [O iii] CELs (Te(0 ii 
Vl/Fll. IPeimbert et al. 2014h . Additionally, for 


A 46 we co uld measure the Tg sensitive ra tio O ii 
4089/4649 (IWesson. Liu, fc B^arlowl [200^ . As it 
can be seen in Table El all these diagnostics give 
values of the electron temperature that are much 
lower than those derived from CELs diagnostics. 
All this clearly points to the presence of a nebular 
component of cold ionized gas where most of the 
ORL emission a rise, as it has been claimed by e.g. 
Liu et al. ( 2006ll . 


Using He+/H+, He^“''/H+, and 0^+ /H+ ratios 
from ORLs (see Table [7]), and 0'''/H+ from ORLs 
in Ou5 and from CELs rescaled to ORLs using 
0'''/0^+ CELs ratio for A 46 and A 63 (see Ta- 
blelT]), the contribution to [O iil]A4363 is estimated 


4.2. Chemical abundances 


Ionic chemical abundances from CELs were 


computed using PyNeb (ILuridiana et al.l 1201411 


and the atomic data in Table ID Errors in the 
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line fluxes and the physical conditions were prop¬ 
agated via Monte Carlo simulations. The average 
rie and Te([0 ill]) for the whole nebula computed 
in the previous section were used. Assuming other 
density values (e.g. the one measured for a spe¬ 
cific ion) would cause very small changes on the 
majority of ions, with the exception of O'*', S+, 
Cl^’*' and Ar^+, whose abundances have a more 
marked dependence on rig. However, even so the 
main conclusions of this work would not change. 
Ionic abundances are presented in Tabled 
As with recombination lines. He’*' and 
abundances were computed from the He i 4471, 
5876, 6678 A and He ii 4684 A lines using the 
updated atomic data in Table [S] In our spec¬ 
tra, we detected and measured several heavy el¬ 
ement ORLs, mainly of O ii and C ii, but also of 
O I. To compute 0^+/H+ ratio we used multi- 
plet 1 O II ORLs around 4650 A, which are the 
brightest and widely use d O ii ORLs (see e. g. 


Garcfa-Roias et al. 20131 and references therein). 
We also computed the C^+/H+ ratio using C ii 
at 4267 A; this is a Sd — 4/ transition that is, in 
principle, purely excited by recombination. The 
sources of the recombination coefhcients that we 
have used to compute the ioni c abundances of C^~* ~ 
and 0^~*~ are listed in T able[5] iTsamis et al. ( 2003 ) 
and Ruiz et al. ( 20031 ) pointed out that the upper 
levels of the transitions of multiplet 1 of O ii are 
not in local thermodynamic equilibrium (LTE) for 
densities ne<10"‘ cm“^, which is the case for our 
three objects, and the abundances derived from 
each individual line may differ by factors as large 
as four. To recover the correct abundances for in¬ 
dividual lines, we applied the n on-LT E corrections 
estimated by Peimbert et al.l (12005 ) to our data. 
We also computed the 0^“''/H+ ratio from the sum 
of all li nes of the mul t iplet following the recipe 
given by Esteban et al.l dlOQSl) . This method takes 
into account the intensity of the whole multiplet 
and is not affected by non-LTE effects. This is the 
value that we finally adopt as representative of the 
02+/H+ ORLs abundance ratio (Table [T]). 


Additionally, we computed the rig that mini¬ 
mizes the dispersion of individual abundances ob¬ 
tained from the O ii multiplet 1 lines. This quan¬ 
tity, labelled as ne(0 ii), is shown in Tabled For 
A 46 and A 63, ne(0 ii) is very similar to ne([0 ilj). 
In Ou5, ne(0 ii) is higher than electron densities 
derived from CELs, in agreement with what found 


in Hf 2-2 by iLiu et al 


(I20n6h . 


Finally, in Ou5 we also computed the 0+/H+ 
ORLs ratio (see Table [T]) , as we detected an emis¬ 
sion line that we identify as the blend of three 
O I ORLs at ~7772 A (see Tab[2|). These lines 
are hardly detected in deep emission spectra of 
other PNe and H ii regions, mainly because they 
are intrinsically faint and frequently blended with 
strong telluric emission lines (which is not the case 
in the spectrum of Ou5). 

The O^^ ad/was computed for each object from 
the comparison of the ORLs and CELs ionic abun¬ 
dances. In Ou5, we could also compute the ad/for 
0+. Results are shown in Table [7| It is clear 
that the adjs of these nebulae- and especially in 
A 46 and Ou5- are very large, independently of 
adopting or not a correction for recombination in 
the calculation of Tg. A detailed discussion on the 
possible origins of the large adjs of these objects, 
the largest ones found in PNe and two orders of 
magnitude larger than in H ii regions, is presented 
in Sect. H) 


4-2.1. Total abundances 

Total abundances for N, C, O, Ne, S, Ar and 

Cl were computed using th e ionization corrections _ 

factor s (ICFs) published bv IPelgado-Inglada. Morisset fc Stasihska 
( 20141 ) from a large grid of photoionization models 
of PNe. The He abundance was computed by sim¬ 
ply adding the He+ and He^"*" abundances. Total 
abundances are presented in Table [S] 


4.3. Spatial variations 

The low surface brightness of these nebulae 
does not allow us to study in detail spatial varia¬ 
tions of the chemical abundances and adjs. Useful 
insights can however be gained by inspecting the 
spatial distribution of selected CELs and ORLs 
along the spectrograph slit. They are presented in 
Figs, m and [S] for A 46 and Ou5. To decrease the 
noise in the faintest emission lines, profiles were 
smoothed with a boxcar average of 5 pixels (1") 
width for A 46, and 3 pixels (0".6) for Ou5. Pro¬ 
files have been arbitrarily scaled in flux to plot all 
lines using the same display. A 63 is not shown 
because of the low S/N of several of these lines. 

In the light of Ha (Figure [Ij, A 46 looks amor¬ 
phous and diffuse. Even so, its brightest regions 
seem to trace the outline of a low surface bright- 
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Fig. 4.— Spatial profiles along the slit of selected 
emission lines for A 46. The cc-axis zero point 
corresponds to the part of the slit closest to the 
central star. The emission peak for the ORLs C ii 
and O II lines is marked in Figure [TJ 



Fig. 5.— Same as in Figure HI but for Ou5. The 
central peak is the continuum emission from the 
central star. The [O ii] 3727 profile is not shown 
here because of the significant, asymmetrical con¬ 
tamination of the background Galactic emission. 


ness inner “cavity” with a bipolar shape whose 
long axis is approximately oriented at position an¬ 
gle P.A.=45°. Figure H] shows that Hi, He i and 
strong CELs such as [O ill] 4959 have a similar 
spatial distribution that follows what is seen in 
Figured] On the other hand, the ORLs emission, 
illustrated by C ii 4267 and O ii 4649, peaks in the 
innermost regions, inside the low surface bright¬ 
ness Ha “cavity”. Their peak is indicated by a 
green cross in Figured! This central enhancement 
is a standard characteristic of the ORLs abun¬ 


dances (I^ 2003 : Tsamis et al. 20081) . Most of 
the fluctuations observed in the ORLs spatial pro¬ 
files are noise in these faint lines, which prevents 
the detection of small-sca le spatial varia t ions s uch 
as for instance done by iTsamis et al. ( 2008f) in 
NGC 7009. More puzzling is the behaviour of the 
auroral line [O ill] 4363, which shows a very similar 
spatial distribution to the ORLs, even if its peak 
is displaced to the west by a couple of arcseconds. 
The marked difference between the [O ill] 4959 and 
[O III] 4363 spatial profiles apparently indicates 
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that III]) decreases from 20 000 K in the in¬ 

nermost regions to 11000 K in the outer parts. 
However, in the previous section we have shown 
that auroral lines can be significantly affected by 
recombination. Table [5] provides some crude es¬ 
timates for the integrated spectrum, but the ac¬ 
tual contribution of recombination to each line as 
a function of position can only be roughly esti¬ 
mated with the present data. The spatial profile of 
[O III] 7330, whose emission should be dominated 
by recombination, is also very similar to that of 
the ORLs and qualitatively supports the hypoth¬ 
esis that at least part of the apparent Te([0 ill]) 
radial gradient could be spuriously produced by 
variable contribution of recombination to the au¬ 
roral lines. 


This is confirmed by the following exercise. We 
extracted the spectrum within 7 arcsec of the 
ORLs peak in the spatial profiles of A 46. Even if 
some line fluxes are one order of magnitude smaller 
than in the integrated spectra, most of the criti¬ 
cal lines could be measured. By proceeding as 
for the integrated spectra, 48% of the [O ill] 4363 
flux is estimated to be produced by recombina¬ 
tion and was corrected accordingly. The value 
Te([0 iii])=17800±1900 K that would be obtained 
by neglecting recombination is then lowered to 
12950±1150 K when it is instead taken into ac¬ 
count. A Te(0 II Vl/Fl) of 3840±200 K is com¬ 
puted, which is only slightly lower than for the 
integrated spectrum. Densities are about twice 
as high, namely ne([S ii])= 560l4gQ° cm“^ and 
ne([0 Ii])=5750t295o cm“^. Total abundances 
determined from CELs only vary slightly com¬ 
pared to those found using the integrated spec¬ 
trum. But, as far as the abundance is con¬ 
cerned, the CELs abundance slightly decreases to 
7.70±0.11 while the ORL one more markedly in¬ 
creases to 10.2ltg;4g. This results in an ac(/(0^+) 
over 300 in the inner regions of A 46! 


Ou5, which has a clean bipolar geometry 
( Corradi et al. 2014ll . shows similar results (Fig¬ 
ure [S]). Hi, He i, [O iii]4959, and [N ii]6583 have 
a similar spatial distribution, peaking at the limb- 
brightened edges of the bipolar lobes. On the 
other hand, ORLs such as C ii and O ii, as well 
the auroral lines of [O ill] and [O il] are more 
concentrated toward the center, and are mutually 
similar. Therefore recombination should be care¬ 
fully considered when the auroral lines are used to 


determine the physico-chemical conditions in this 
kind of nebulae. 


5. Total Ha flux and nebular mass 


A 63 and Ou5 were observed by the IPHAS 
Ha photometric survey of the Galactic plane 
( Drew et al.ll2005 1. Therefore their IPHAS Ha-|-[N ii] 
images can be flux calibrated using the in- 
formation in the IPH AS Second Data Release 
(jBarentsen et al.l 120141 1. The observed, total Ha 
fluxes of these nebulae, after removing the contri¬ 
bution of [N il] 6548,6583 to the IPHAS Ha filter 
using our spectroscopy, are listed in Table [1] Er¬ 
rors are of the order of 10%. Taking advantage of 
the fact that the image of A 46 in Figure [T] was 
obtained with the same instrumentation and filter 
in a clear night, we adopted the mean value of the 
zeropoint of the IPHAS survey to also calibrate 
this image. The additional error introduced by 
this assumption is another 10%. 


With the total Ha flux and the physical param¬ 
eters determined from spectroscopy, and adopt¬ 
ing the distances quoted in Table [U the nebular 
ionized mass can be estimated. The total ionized 
mass of hydrogen was computed by applying the 
expression: 




471 7^(LI/3) nip 
/uzh/3 ne {H°,Te) 


where F(H/3) is the dereddened H/3 flux, D is the 
distance to the object, nip is the proton mass, rie 
is the electron density, is the energy of an 

H/3 photon and, finally, (i7°, Te) is the effec¬ 
tive recombination coefficient of H/3. The usual 
dependence on the square of the distance appears 
in the formula. In this respect, it should be noted 
that the distances of A 46 and A 63 were origi¬ 
nally computed through the stellar parameters of 
the central binaries, which should provide more re¬ 
liable determinations than any statistical method, 
especially when the peculiar nature of these neb¬ 
ulae is considered. 

As (xTe the total mass is 

roughly proportional to the adopted electron tem¬ 
perature Te and inversely proportional to the elec¬ 
tron density rig. The dual-component nature of 
these nebulae prevents a precise mass determina¬ 
tion, as the relative amounts of H i emitting at 
the low Te indicated by the Balmer jump and the 












He I or O ii ORLs, and that emitting at the much 
higher CELs Te is not known. Strictly speaking, 
at least for A 46 one would be tempted to consis¬ 
tently adopt for hydrogen the temperature com¬ 
puted using the Balmer jump, namely 1000 K. As 
with densities, there is also a clear difference (by 
up to a factor of 8) between those computed from 
[O ii] and those computed from [S il]. Assum¬ 
ing Te=1000 K and the average rie from CELs, we 
determine an H+ mass of 1x10“^, 5x10“^, and 
5x10“^ Mq for A 46, A 63, and Ou5, respectively. 


However, even if the total amount of metals in 
the hot and cold gas phases is c omparable as fo r 
instance estimated for Hf 2-2 ( Liu et al. 2006[l . 
the much higher metallicity of the cold compo¬ 
nent implies that the majority of hydrogen must 
be in the hot component. It seems therefore rea¬ 
sonable to adopt re([0 III]), which would imply 
H+ masses almost one order of magnitude larger, 
namely 8x10“^, 4x10“^, and 5x10“^ Mq for 
A 46, A 63, and Ou5, respectively. 


Einally, in the less likely case that all hydro¬ 
gen comes from a plasma with re([0 ill]) and the 
lowest electron density measured (ne([S ii])), we 
obtain 4x10“^, 1x10“^, and 2x10“^ Mq, re¬ 
spectively. These should be just taken as con¬ 
servative upper limits to the total amount of 
H+ in these neb ulae (for the adopted distances). 


Liu et al.l (120061) found a similarly low H+ mass of 
1.4x10-2 Mq for Hf 2-2. 

Therefore, in spite of the large uncertainties 
in these calculations, the clear conclusion is that 
these nebulae with large adjs and close binary 
central stars have a very small total ionized mass, 
one or more orders of magnitude s lower than typ - 
ical of Galactic PNe (see e.g. IPottascb 19841) . 
The total nebular mass could in principle be 
larger if part of it is in neutral/molecular form. 
However, spectral indic ators such as the absence 
of [O i] 6300 emission ( Harman fc Seatoiil Il966t 
Kaler fc Jacobv 19891) . which is thought to be 


produced at the interface between ionized and 
neutral gas via charge exchange reactions, point 
to optically thin, density bounded nebulae. 


6. Discussion 


PNe. In particular, A 46 has the largest ad/ever 
measured, with the exception of the hydrogen- 
defi cient knots that are directly o bserved in Abell 
30 ( Wesson. Liu, fc Barlow! [200 . Our third tar¬ 
get, A 63, has a lower adf (~10), but still above 
the average value in PNe. 

Our three target PNe have close binary cen¬ 
tral stars with orbital periods around ten hours. 
Therefore these observa tions, together with the 
additional case of Hf 2-2 ( Liu et al]l2006 ). reveal a 
clear connection between the most extreme repre¬ 
sentatives of the abundance discrepancy problem 
and the binarity of the central sources. The above- 
mentioned Abell 30 (and the morphologically sim¬ 
ilar Abell 58 aka V605 Aquilae, which experienced 
a nova-like outburst in 1 905) were also prop osed 
to contain a close binary ( Wesson et al.ll2008l) . but 
no direct evidence has been found yet. 

Our spectroscopic analysis supports the previ¬ 
ous interpretation that two diff erent gas phases 


coexist in these nebulae (e.g. Liu et ahl 2006t 


Tsamis et al.ll20(38l) : hot gas at 10'* K with stan¬ 


dard metallicity where the CELs can be efficiently 
excited, and a much cooler (10^ K) plasma with a 
highly enhanced content of heavy element abun¬ 
dances (which cause the cooling) where only ORLs 
form. How much each component contributes to 
the total mass, and how they are distributed and 
mixed, is poorly known, because existing data are 
not able to spatially resolve them. One common 
characteristic, also shared by the nebulae stud¬ 
ied in this work, is that the ORLs emission is 
enhanced in the central regions of the nebulae, 
where the adf r eaches a value over 300 in A 46. 


Liu et al] (200^ estimated that similar amounts 
of metals are in the hot and cold gas of Hf 2-2, 
but the amount of H and He in each component 
cannot be determined. This prevents the calcula¬ 
tion of precise X^~*~ /H~*~ abundanc es for each gas 
phase individually (jLiu et al.ll2006l) , and therefore 
all values in Table |S] should be taken with caution. 
Eor instance, the large helium content that is mea¬ 
sured in A 63 and Ou5 may not be com pletely real 
(see also Liu 2003t iTsamis et al.ll2008h . Only ionic 
abundance ratios, such as for example N+/0+, are 
free of this bias. In this respect A 46, and to a 
lesser extent A 63 (because of its large errors) and 
OufH, have low N+/0+RiN/0 CELs abundance 


Two of the three PNe that we studied, namely 
A 46 and Ou5, have abundance discrepancy fac¬ 
tors that are among the largest ones found in 
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^The much lower N+/0+ ratio for this nebula computed by 



























































ratios, typical of low-nietallicity SMC PNe, or al¬ 
ternately indicating substantial oxygen enhance¬ 
ment compared to nitrogen. But even this ratio 
should be taken with caution, considering the high 
excitation of the nebulae. We have also noticed 
that CELs temperature gradients within a nebula 
can be at least partly spurious, caused by a vary¬ 
ing contribution of recombination to the emission 
of the auroral lines. This additional uncertainty 
also affects the determination of the CELs abun¬ 
dances. 

6.1. The nature of these nebulae and the 
implications for the abundance dis¬ 
crepancy problem 

The main result of this work is that any ex¬ 
planation of the abundance discrepancy problem 
must include the fact that the most pathological 
cases are planetary nebulae with a close binary 
central star. This casts light on the origin of the 
problem. Let’s consider previously proposed and 
new scenarios, and their pro et contra. 

A key point to understand the evolution of 
these systems is the small ionized mass of the neb¬ 
ulae (Sect. [5]). Let’s take 10“^ Mq as a repre¬ 
sentative figure of the nebular mass for the sys¬ 
tems that we have considered. The total mass 
in the high-metallicity component would then be 
a (likely small) fraction of it. These low masses 
dispute the usual interpretation of the origin of 
these nebulae. It is generally thought that PNe 
with close binary central stars went through a 
common-envelope phase when the radius of the 
primary star increased in the AGB phase and en¬ 
gulfed the companion. The CE interaction caused 
the shrinkage of the orbit to the current small 
separations (~2 R 0 for the four objects in Ta¬ 
ble [T|) and provoked the complete ejection of the 
envelope, forming the PN. However, the nebular 
masses that we have estimated are far too low 
compared to the AGB envelope masses expected 
for 0.5-0.6 M 0 cores and standard initial-to-final 
mass relationships, e ven if some mass loss p rior to 


the AGB is allowed ( 

De Marco et al.l 2011 

. The 

exercise has been done by Afsar & Ibanoglu 

(2nnsh 

and De Marco et al. ( 

2011) for A 46 and A 63. for 


ICorradi et al.l <120141) was caused by the use of the auroral 
[O ll]7320,7730 lines, that we have shown here to be severely 
affected by recombination processes. 


which stellar pa r amete rs are known (Table [T]). In 
De Marco et al. ( 201ll) . consistency between the 


stellar evolution and binary parameters is reached 
for CE masses of 0.4 M 0 for A 46 and 1.2 M 0 for 
A 63, i.e. one to three orders of magnitude larger, 
respectively, than what we measured. The similar 
or only slightly larger nebular masses of Ou5 and 
Hf 2-2 add strength to this conclusion, even if the 
stellar parameters for these two systems are not 
precisely knowr0. This missing-mass problem is 
even more severe considering that the whole CE 
is thought to be eje cted in a very short tim escale 
(^few years, see e.g. lSandauist et aPIlQQSll . limit¬ 
ing the possibility that the PN is dispersed over a 
large volume. 


As it is evident that the systems at some point 
went through a CE phase (given the present small 
orbital separation), the envelope of the primary 
star would have been ejected a long time ago 
and would have now vanished in the interstellar 
medium. Note that the 0.5 M 0 mass of the cen¬ 
tral star of A 46 is at the limit between RGB and 
AGB core masses. It cannot be excluded that in 
this case (but can be ruled out for A 63 which has 
a typical AGB core mass of 0.63 M 0 ) the CE in¬ 
teraction, ejection, and orbital shrinking occurred 
in the RGB phase of the primary. 


Then the currently observed low-mass nebu¬ 
lae would be produced much after the CE phase, 
in post-CE mass loss episodes when the massive 
AGB envelope has already dispersed. Even so, 
their central stars have to be still hot and lumi¬ 
nous to be able to ionize the nebulae, as it is indeed 
observed (see e.g. Afsar fc Ibanoelu 2008lj . 


6.1.1. Very late thermal pulses (VLTPs) 


This hypothesis has been discussed in se veral 
articles ( Liu et all 2006t Wesson et al. 2008 ). In 
general, thermal pulses in the post- AGB phase ar e 
expected to be relatively frequent (lBlockeiil200ll) . 
In the best studied case, Sakurai’s o bject, the esti¬ 
mate d ejected mass is 2x10“^ M 0 (IHaiduk et al. 
20051) . i.e., in the possible range of what we have 


measured. Gas falling back from the ejected enve¬ 
lope, accreting on the post-AGB star, and perhaps 
trig gering a VLTP was pr e dicted and discussed by, 
e.g., Frankowski fc Soker ( 20091) . In those specific 
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®For Ou5, they are very likely in the range of A 46 and A 63 
(Corradi et al. 2015, in preparation). 



























































calculations, however, the total amount of fallback 
material, 10“^ to 10“^ M©, seems on the low side 
of what we observed, and the timescale for accre¬ 
tion is too short to allow the massive, ejected CE 
to vanish completely. 

The VLTP ejecta are expected to basically rep¬ 
resent the intershell abundances, with C/0^1. In 
the case of Ou5, the ORL C/0 ratio could be esti¬ 
mated as C/0=a where a=1.12 is the 

ICF f rom [Delgado-Inglada. Morisset &: Stasihska 
(1201411 that could be consistently computed using 
the ORL 0+ and 0^+ abundances. We obtain 
0/0=0.55±0.24 for this nebula. For A 46 and 
A 63, no ORL O'*" abundance determination is 
available, and therefore the ICF can only be es¬ 
timated using CELs which would indicate values 
very close to unity. We obtain C/0=0.78±0.11 
for A 46, and 1.66±0.47 for A 63 which is the only 
nebula where the C/0 ratio seems to be above 
unity. Therefore, in spite of the uncertainty in¬ 
volved in the use of the faint ORLs, the low C/0 
ratios e stimated in A 4 6 and Ou5, as well as in 
Hf 2-2 (iLiu et al.ll200il . argue against VLTPs as 
the origin of the high-metallicity ejecta causing the 
large observed adjs. 

6.1.2. Nova-like outbursts 

Nova outbursts are natural phenomena in in¬ 
teracting binaries. They result from thermonu¬ 
clear runaways on the surfaces of white dwarfs 
(WDs) accreting hydrogen-rich ma tter in binary 
systems (see e.g. ICehrz et al.lll998h . Nova ejecta 
are metal enriched by mixing between the ac¬ 
creted envelope and the WD core. Therefore, in 
principle, the metal-rich gas in our target neb¬ 
ulae could be ejected in nova outbursts from 
the close binary centr a l star s, as suggested by 
Wesson. Liu, fc Barlowl ( 2003h and Wesson et al 


( 2008l l for A 30 and A 58. In these two cases, the 


hypothesis of a nova eruption could explain the 
large neon enrichment of their hydrogen-deficient 
knots, if produced in a massive (M>1.2 Mq) 
ONeMg WDs. For A 58, this scenario (in com¬ 
bination with a final helium shell flash) or the 
alternati ve possibility of a stell a r mer ger are eval¬ 
uated in Lau. De Marco fc Liu ( 20Illl . 

However, there is no indication that mass trans¬ 
fer from the secondaries to the primaries is cur¬ 
rently going on in the systems that we have stud¬ 
ied, because the secondary stars - albeit inflated 


- do not hll their Roche lobe (jAfsar &: Ibanoglu 
20081 1. To feed the WD, accretion must therefore 
come from circumbinary material left-over in the 
systems from the CE phase. A possibility is that 
the envelope is not removed completely in the CE 
interaction and part of it remains gravitationally 
bounded to the system. The non-accreted part 
of this circumstellar gas, with normal AGB abun¬ 
dances, could then interact with the nova ejecta 
and partially mix with it producing the observed 
dual-component nebulae. In this respect, it should 
be noted that the existing simulations of the CE 
process fail to eject the entire envelope, a large 
fraction of wh ich remains grayitationallY bound 
to the system (ISandauist et al. I998t Passv et "ni 
2012 ; Ricker &: TammI 2012 ). This may result in 
a series of CE and fail-back events, and in the 
forma tion of a left over, l o w mass circumb inary 


disk ( Kashi fc Sokei 2011 : iKuruwit^ l2015l Ku 


ruwita, De Marco & Staff, in preparation), which 
could in turn form metal-rich dust as perhaps ob- 
served in the putativ e post-CE system V1309 Sco 
( Nicholls et al.l2013ll . 

Possible problems are the very low mass of clas¬ 
sical nova ejecta (<10~^ Mf?^, but generally much 
smaller, see e.g. ICehrz et al. 1998ll . and the short 
timescale (a few years) of their nuclear burning 
and high luminosity phase that are needed to light 
up the PNe. Long-lasting outbursts at high WD 
luminosity can be achieved if the mass accretion 
rate is higher than in classical novae, as for in¬ 
stance occurs in symbiotic (very slow) novae. The 
relevant parameters space should be explored in 
detail to see how this hypothesis ca n fit our post- 
CE PNe (cf. e.g. Yaron et al. 2005). 

An additional important parameter is the ex¬ 
pansion velocity of the ejecta. In the PNe that we 
have studied, the ORL emitting gas phase does not 
show the large expansion velocities (few to many 
hundred km s“^) typical of classical and symbi¬ 
otic novae, which would be easily revealed at the 
resolution of our spectra. While more precise in¬ 
formation is not available yet, a further constraint 
for the modeling is that both the hot and cold gas 
components have expansion velocities more typi¬ 
cal of PNe than of novae. 

6.1.3. Planetary material 

A source of the metal-rich component in the 
nebulae might be planetary debris that survived 
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the CE phase. Indeed, about a quarter of single 
WDs have considerable amounts of metals in their 
photosphere, that are explained by accretion of 
tidally disrupted planetesimals/asteroids/comets 
that survived through all previous evolutio nary 
phases of the star (see e.g. lFarihi et al.ll2014ll . 

The possible contribution of solid bodies to 
the metal- rich c ompo nent in PNe was first men- 
tioned bv iLiul (1200311 and then investigated by 
Hennev fc Stasihska ( 2010ll . They conclude that 
solid bodies could provide the source of the high- 
metallicity gas only if they are meter-sized or 
larger, their evaporation process extends back to 
the final AGB evolution, mixing with the gas com¬ 
ponent at “normal” AGB metallicity is inefficient, 
and systems possess a much more massive popula¬ 
tion of comets than is found in our Solar system. 

Note also that Tsamis et al.l ( 2011 . 201, Sh and 
Mesa-Delgado et al. ( 2012ll proposed that photo¬ 
evaporating protoplanetary disks are the cause of 
the smaller adjs found in H ii regions. 


6.1.4- Planets destruction 

Another possibility is that these post-CE nebu¬ 
lae are produced by the destruction of one or more 
circumbinary Jovian planets. There is increasing 


evidence that post-CE planets exis t jlParsoi^.et^ 


Bear fc Soker 


2014 : IZorotovic fc Schreibeil 201,11 : 

201411 . It is therefore possible that planets in 


unstable orbits spiral in, are tida lly disrupted, 
and c ollide with the post-AGB star ( Bear fc Soked 
20121 1. The hydrogen-rich envelope of the planet 
would be stripped first, and then its metallic core. 
Part of this material would be accreted by the 
WD: at high accretion rates, the gas cannot be 
processed and inflates a red-giant envelope, pro¬ 
ducing a new CE phase and eventually the ejec¬ 
tion of a new nebula with the dual nature (H-rich 
from the planet envelope and metal-rich from its 
core) that we have observed. The new CE would 
also favour the formation of highly bipolar neb¬ 
ulae such as Ou5 and A 63, by forcing the CE 
ejection toward the orbital plane, and/or produc¬ 
ing jets from one of the two stars (as observed 
in A 63). A possible problem with this hypothe¬ 
sis is the apparently enhanced helium abundance 
of our target PNe (Table [8]) , which is not charac¬ 
teristic of planets. However, the He abundances 
in thes e high-adf objects y e likely to be overesti¬ 
mated ( Tsamis et al. 2008h . 


While this scenario should be taken as highly 
speculative at this stage, it has several positive as¬ 
pects such as the fact that the process for planet 
collision is long enough (>10® year) to allow the 
CE to vanish, the ejection could be split into dif¬ 
ferent phases (the planet’s atmosphere and core), 
and the observed nebular masses are in a plausi¬ 
ble range corresponding to 1-10 Jupiter masses. In 
such case. Sir William Herschel would have been 
right in coining the name planetary nebulae! 


6.2. Other implications and conclusions 


If our hypothesis of an intimate relation be¬ 
tween binarity and large adjs is correct, we pre¬ 
dict that many - if not all - the PNe where large 
adjs have been measured are post-CE binaries. 
One already confirmed case is NGC 6778 which 
i ndeed contains a cl ose binary with Por6=0.15 
( Miszalski et al.l[201l[) . Other primary candidates 
to search for binarity are M 1-42, M 2 -36, M 3- 
26, M 3-32. NGC 61,13 and NGC 1501 (Ijj^l2n03 : 
Tsamis et al.l 120041 : lErcolano et al. 20n4ll . as well 


as the already mentioned Abell 30 and Abell 58. 
Note that our conclusion does not imply the op¬ 
posite hypothesis, i.e. that all post-CE PNe have 
large adjs. For exa mple, NGC 5189 has a close bi¬ 


nary central star (iManick. Miszalski fc McBride 
l201,5l l but a low adf ( Carcfa-Roias et al.1 201,lll . 
Also, no oxygen or carbon recombination lines are 
detected in relatively deep spect ra of the Necklace 
binary PN ( Corradi et al. 2011 ). 

If some of the post-CE nebulae are not gen¬ 
uine AGB envelopes, but rather post-CE mass loss 
events or post-RGB systems, the total fraction of 
close binary PN central stars w ould be smaller 
than presently estimated (^15% iMiszalski et al 


2009f) . and perhaps more consona nt with the frac¬ 


tion o f main-sequence bin aries (jPe Marcoetjd 


2013; iDonchin et al.l 12014 but also see [Boffin 


20141) . 


Concluding, we have added new ingredients to 
the discussion of the nature of PNe with close 
binary central stars, the CE evolution, and the 
abundance discrepancy problem. The main result 
of this work is the clear link between binary stellar 
evolution and the formation of high adjs in PNe. 
At the moment, none of the proposed scenarios 
to explain this result clearly stands out, but this 
is likely because observations provide still limited 
constraints. To further progress, a larger sample of 
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objects should be analyzed to highlight the basic 
relationships among the numerous parameters in¬ 
volved. In particular, more precise determinations 
of crucial quantities such as the nebular masses 
and chemical compositions for each (normal and 
enriched) nebular components, and the binary pa¬ 
rameters (to better reconstruct the evolution of 
the systems), are needed. 

Finally, even if the binary nature of the central 
stars seems to be an essential ingredient to under¬ 
stand the abundance discrepancy problem in PNe, 
it is hard to envisage its relevance to other astro- 
physical contexts such as H ii regions. This may 
imply that different causes should be sought to 
fully understand this lo ng-standing astrophysical 
proble m, as proposed bv iGarda-Roias fc Esteban 
( 20071) . On the other hand, as discussed in the 
previous sections, some related phenomena may 
be involved, like for instance the possible contri¬ 
bution of protoplanetary disks in H ii regions and 
of planetary debris in PNe. 

This work is based on observations obtained 
with the 4.2m WHT and 2.5m INT telescopes of 
the Isaac Newton Group of Telescopes, operating 
on the island of La Palma at the Spanish Obser¬ 
vatories of the Roque de Los Muchachos of the 
Institute de Astrofisica de Canarias. The WHT 
spectra of 14 October 2014 were obtained in ser¬ 
vice time. Also based on data obtained from the 
ESO Science Archive Facility. We are extremely 
grateful to Noam Soker and Orsola De Marco for 
their significant contribution, in terms of origi¬ 
nal ideas and proper criticism, to the discussion 
in Sect. O We also thank Valentina Luridi- 
ana for her help with PyNeb, Grazyna Stasihska, 
Christophe Morisset, and Gesar Esteban for a 
critical reading of the manuscript, and the referee 
for his/her very useful comments. This research 
has been supported by the Spanish Ministry of 
Economy and Competitiveness (MINECO) un¬ 
der grants AYA2012-35330, AYA2011-22614, and 
AYA2012-38700. JGR acknowledges support from 
Severo Ochoa excellence program (SEV-2011- 
0187) postdoctoral fellowship. PRG is supported 
by a Ramon y Cajal fellowship (RYC2010-05762). 
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Table 1 

Basic properties of the PNe, central stars, and observational details. 


Name 

PN G 

R.A. (J2000) Dec 

Slit P.A. 
(deg) 

D 

(kpc) 

^ orb 

(days) 

rrics 

(Mq) 

m2 

(Mq) 

i 

n 

F(Ha) 

(erg cm~^ s~^) 

(Mq) 

Abell 46 

055.4+16.0 

18 31 18.29 

+26 56 12.9 

100 

1.7 

0.47 

0.51 

0.15 

80 

8.87x10-12 

10-3 to 10-2 

Abell 63 

053.8-03.0 

19 42 10.20 

+17 05 14.4 

300 

2.4 

0.46 

0.63 

0.29 

87 

2.29x10-12 

10-1 to 10-2 

Ou5 

086.9-03.4 

21 14 20.03 

+43 41 36.0 

297 

5 

0.36 

m2/m.cs 

RiO.4? 

90 

0.91x10-12 

10-3 to 10-1 

Hf 2-2 

005.1-08.9 

18 32 30.93 

-28 43 20.5 


4.25 

0.40 



25-55 


RilO-2 


Note. —Fo r the original sources of the orb ital periods, se e http: / /www.drdiones.n et/?q=node/6. Di stances, binary and stellar parameters for 
A 46 are from |BelL PoIIacco Hilditchl lll994l~) , for A 63 from lAfsar Sc Ibanoglul (120081) . and for Ou5 from ICorradi et al.fll2014l , 2015 in preparation). 
Hf 2-2 is no t directly analize d in this paper, but it is shown here becaus e is the other key target for discuss ion in Sect.[6| Its spectroscopic analysis 
was done bv iLiu et al.l (l2006lV the quoted central star inclination is from l^haub. Bodman &; HillwiS ll2012t) . 

^For the adopted distance (col. 6), and computed as described in Sect. [5] 


























Table 2 

Line identifications, observed fluxes and their % errors 


Ao (A ) 

Ion 

Mult. 

Abell 46 

Abell 63 


Ou5 

F(A) 

Error (%) 

F(A) 

Error (%) 

F{A) 

Error (%) 

3634.25 

He I 

28 

1.21 

16 





3679.36 

Hi 

H21 

0.80 

23 





3682.81 

Hi 

H20 

0.67 

26 





3686.83 

Hi 

H19 

1.03 

18 





3691.56 

Hi 

H18 

1.51 

14 





3697.15 

Hi 

H17 

0.95 

20 





3703.86 

Hi 

H16 

2.56 

10 





3711.97 

Hi 

H15 

2.49 

10 





3721.83 

[S HI] 

2F 

2.02 

11 



2.42 

31 

3721.93 

Hi 

H14 

* 

* 



* 

* 

3726.03 

[O II] 

IF 

23.91 

6 

11.01 

11 

26.01 

8 

3728.82 

[O II] 

IF 

14.82 

6 

6.61 

12 

20.68 

8 

3734.37 

Hi 

H13 

2.20 

10 

2.97 

18 

3.16 

25 

3750.15 

Hi 

H12 

3.88 

8 

3.78 

16 

3.12 

25 

3770.63 

Hi 

Hll 

4.35 

7 

3.64 

16 

4.17 

19 

3797.63 

JS HI] 

2F 

5.49 

7 

6.56 

12 

5.16 

16 

3797.90 

Hi 

HIO 

* 

* 

* 

* 

* 

* 

3819.61 

He I 

22 

1.39 

14 





3835.39 

Hi 

H9 

6.60 

6 

7.15 

11 

7.26 

12 

3868.75 

]Ne III] 

IF 

39.43 

5 

25.37 

9 

74.32 

6 

3888.65 

He I 

2 

22.66 

5 

22.78 

9 

22.17 

7 

3889.05 

Hi 

H8 

* 

* 

* 

* 

* 

* 

3964.73 

He I 

5 

0.91 

19 





3967.46 

[Ne III] 

IF 

12.61 

5 

7.33 

10 

21.00 

7 

3970.07 

Hi 

H7 

16.69 

5 

17.41 

8 

16.57 

7 

4009.26 

He I 

55 

0.67 

25 





4026.21 

He I 

18 

3.58 

7 

2.86 

16 

3.59 

19 

4068.60 

JS II] 

IF 

4.38 

6 



4.30 

16 

4069.62 

O II 

10 



1.69 

24 



4069.89 

O II 

10 



* 

* 



4076.35 

JS II] 

IF 

2.56 

8 



2.99 

21 

4078.84 

O II 

10 

0.37 






4085.11 

O II 

10 

0.67 

25 





4089.29 

O II 

48 

1.16 

15 





4097.22 

O II 

20 

2.69 

8 



3.20 

20 

4101.74 

Hi 

H6 

26.32 

4 

25.87 

7 

26.43 

6 

4104.99 

O II 

20 

0.65 

25 





4107.09 

O II 

62 

0.50 

32 





4119.22 

O II 

20 

1.20 

15 





4132.80 

O II 

19 

0.65 

25 





4153.30 

O II 

19 

1.12 

16 





4267.15 

C II 

6 

6.45 

4 

3.99 

11 

5.11 

12 

4275.55 

O II 

67 

1.12 

15 



1.52 

35 

4276.75 

O II 

67 

0.99 

17 





4294.92 

O II 

54 

0.56 

28 





4303.61 

O II 

65 

1.57 

11 

0.93 

37 

1.60 

32 

4303.82 

O II 

53 

* 

* 

* 

* 

* 

* 

4317.14 

O II 

2 

0.56 

28 





4319.63 

O II 

2 





1.04 



18 



























Table 2 —Continued 


^0 (A ) 

Ion 

Mult. 

Abell 46 

Abell 63 

Ou5 


F(A) 

Error (%) 

F(A) 

Error (%) 

E(A) Error (%) 

4340.47 

Hi 

H5 

47.34 

3 

48.15 

5 

49.45 

4 

4345.56 

O II 

2 

0.46 

33 



1.10 


4347.41 

O II 

2D-2Do 





1.22 

40 

4349.43 

O II 

2 

0.93 

18 



1.44 

35 

4363.21 

[O III] 

2F 

4.88 

5 

1.07 

32 

6.33 

9 

4366.89 

O II 

2 

0.90 

18 





4379.55 

Ne II 

60 





1.84 

27 

4387.93 

He I 

51 

0.71 

22 





4414.90 

O II 

5 

0.35 






4416.97 

O II 

5 

0.44 

34 





4471.47 

He I 

14 

6.49 

4 

6.43 

7 

6.27 

9 

4491.23 

O II 

86a 

0.80 

20 





4541.59 

He II 

4.9 

0.76 

21 



2.56 

18 

4609.44 

O II 

92a 

0.91 

17 



0.91 


4638.86 

O II 

1 

1.24 

13 

0.59 


1.51 

28 

4640.64 

N III 

2 





1.30 

32 

4641.81 

O II 

1 

2.19 

8 

0.94 

32 

2.68 

17 

4641.85 

N III 

2 





* 

* 

4649.13 

O II 

1 

3.14 

6 

0.75 


3.81 

12 

4650.84 

O II 

1 

1.81 

9 

0.53 


1.80 

24 

4661.63 

O II 

1 

1.51 

11 



1.89 

23 

4676.24 

O II 

1 

0.72 

21 





4685.68 

He II 

3.4 

26.41 

1 

6.83 

6 

71.11 

1 

4711.37 

[Ar iv] 

IF 

1.25 

13 

1.19 

26 

3.61 

12 

4740.17 

[Ar iv] 

IF 

0.93 

17 



2.73 

16 

4859.32 

He II 

4.8 

1.97 

8 



3.75 

11 

4861.33 

Hi 

H4 

100.00 

0 

100.00 

1 

100.00 

1 

4921.93 

He I 

48 

1.58 

10 

1.78 

17 

1.46 

25 

4958.91 

[O III] 

IF 

119.93 

1 

85.47 

1 

184.58 

1 

5006.84 

[O III] 

IF 

361.10 

1 

259.49 

1 

567.04 

1 

5411.52 

He II 

4.7 





5.15 

7 

5517.71 

JCl III] 

IF 





0.70 

37 

5666.64 

N II 

3 





0.54 


5679.56 

N II 

3 

0.66 

21 

0.51 


0.85 

29 

5754.64 

[N II] 

3F 

0.40 

33 

0.47 


0.85 

29 

5875.64 

He I 

11 

22.06 

4 

22.47 

8 

20.12 

6 

6312.10 

[S III] 

3F 





1.00 

21 

6402.25 

Ne I 

1 





0.44 


6406.30 

He II 

5.15 





* 

* 

6461.95 

C II 

17.04 

0.51 

25 



0.54 

34 

6527.11 

He II 

5.14 





0.45 

39 

6548.03 

[N II] 

IF 

1.40 

12 

2.26 

15 

5.44 

9 

6562.82 

Hi 

H3 

318.63 

7 

347.18 

11 

323.61 

8 

6578.05 

C II 

2 

0.63 

21 

0.40 


0.64 

29 

6583.41 

[N II] 

IF 

4.95 

7 

5.26 

12 

15.05 

8 

6678.15 

He I 

46 

6.27 

7 

6.12 

13 

5.72 

9 

6683.20 

He II 

5.13 





0.35 


6716.47 

[S II] 

2F 

1.42 

12 

0.60 

32 

4.32 

10 

6730.85 

[S II] 

2F 

1.28 

13 

0.65 

30 

3.45 

10 
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Table 2 —Continued 


Ao (A ) 

Ion 

Mult. 

Abell 46 

Abell 63 

Ou5 


F(A) 

Error (%) 

F(A) 

Error (%) 

F(A) Error (%) 

6890.88 

He II 

5.12 





0.35 


7005.67 

[Ar v] 

3P-1D 





0.24 


7065.28 

He I 

10 

2.06 

10 

2.30 

16 

1.79 

14 

7135.78 

[Ar III] 

IF 

5.58 

9 

8.33 

14 

13.12 

10 

7177.50 

He II 

5.11 

0.24 




0.71 

23 

7231.34 

C II 

3 





0.35 


7236.42 

C II 

3 





0.74 

22 

7237.17 

C II 

3 





* 

* 

7281.35 

He I 

45 





0.28 


7318.92 

[O II] 

2F 



1.58 

18 

3.26 

12 

7319.99 

]o II] 

2F 



* 

* 

* 

* 

7329.66 

]o II] 

2F 

2.74 

10 

1.24 

20 

2.73 

12 

7330.73 

]o II] 

2F 

* 

* 

* 

* 

* 

* 

7530.54 

[Cl iv] 

IF 





0.32 


7530.57 

C II 

16.08 





* 

* 

7592.74 

He II 

5.10 





1.04 

17 

7751.10 

[Ar III] 

2F 





2.76 

13 

7771.93 

O I 

1 





1.00 

18 

7774.17 

O I 

1 





* 

* 

7775.39 

O I 

1 





* 

* 

8045.63 

]C1 iv] 

IF 





0.74 

21 

8196.48 

C III 

43 





0.37 

33 

8236.77 

He II 

5.9 

0.72 

19 



1.81 

15 

8467.25 

Hi 

P17 





0.54 

24 

8502.48 

Hi 

P16 





0.39 

30 

8545.38 

Hi 

P15 





0.48 

26 

8598.39 

Hi 

P14 



0.77 

26 

0.65 

22 

8665.02 

Hi 

P13 





1.04 

18 

8703.25 

N I 

1 



0.42 

37 



8750.47 

Hi 

P12 

1.02 

16 

1.04 

24 

1.13 

17 

8862.79 

Hi 

Pll 





1.44 

16 

9014.91 

Hi 

PIO 

1.78 

13 

1.69 

22 

1.75 

16 

9068.60 

[S III] 

IF 

8.46 

12 

14.67 

20 

20.74 

15 

9229.01 

Hi 

P9 

1.97 

13 

1.97 

13 

3.06 

15 


Note. —indicates that this emission line is blended with nearby ones. The total flux of 
a blend is indicated in the first listed line. indicates uncertainties larger than 40%. 
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Table 3 

Extinction, temperatures and densities 



Abell 46 

Abell 63 

Ou5 

c(H/3)’^ 

0.22±0.08 

0.55±0.14 

0.94±0.10 

c(H/3)'^ 

U.U^_Q Q4 

0.38±0.11 

0.77±0.14 

Te ([0 III]) 

12750±200 

8750±950 

11900±500 

Te ([0 III])^ 

12050±200 

7400±550 

10150±300 

Te ([N II]) 

>37000 

>24000 

20000±8500 

Te ([N II])" 

>30000 

>22000 

18800±5000 

Te (He I 5876/4471) 

1950ti°5° 

150U_2^qqq 

2800 lt?™ 

Te (He I 6678/4471) 

2150ti«™ 

oocn+SOOO 

ZZDU_i250 

32001^™ 

Te (Baimer jump) 

1150±550 



Te (0 II 4089/4649) 

800 ; 280 o 


<1000: 

Te(0 II Vl/Fl) 

4300±100 

4525±125 

5020±360 

Re ([S Ilj) 

O^U_220 

600 ; 2 ioo 

iOU_ioo 

Re([0 II]) 

071^0+880 

OU_ggQ 

2600™ 

1200;5™ 

He ([Ar iv]) 

<220: 


300t«°° 

He (CELs, adopted) 

1590±600 

iODU_2ioo 

ODU_3qq 

Re (0 ll) 

2960 

1940 

3900 


Note.— Temperatures are in units of Kelvin, and densities 
in cm“^. indicates uncertain values. 

'^Assuming re=12500 K and ne=1000 cm“^. 

'^Assuming Te=1000 K and ne=1000 cm“^. 

After correction for the estimated contribution of recombi¬ 
nation to the auroral lines (see text and Table E]) . 
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Table 4 

Atomic dataset used for collisionally excited lines 


Ion 

Transition probabilities 

Collisional strengths 

N+ 

Froese Fischer & Tachiev ('2004P 

Taval 120111 

0+ 

Froese Fischer & Tachiev 120041 

Kisielius et al. 120091 

0^+ 

Froese Fischer & Tachiev 12004P 

Storev et al. 120141 

Ne2+ 

Galavis et al. ('19971 

McLaughlin & Bell 120001 

S+ 

Podobedova et al. (20091 

Taval & Zatsarinnv 120101 

S2+ 

Podobedova et al. (20091 

Taval & Guota 119991 

CP+ 

Mendoza 119831 

Butler & Zeiooen 119891 

Cl3+ 

Kaufman & Sugar 119861 
Mendoza & Zeinoen 11982al 

Ellis & Martinson 119841 

Galavfs et al. 119951 

Ar2+ 

Mendoza 119831 

Kaufman & Sugar 119861 

Galavis et al. 119951 

Ar3+ 

Mendoza & Zeiooen 11982bl 

Ramsbottom & Bell 119971 

Ar'‘+ 

Mendoza & Zeinoen 11982al 
Kaufman & Sugar 119861 

LaJohn & Luke 119931 

Galavis et al. 119951 


Adopting the A-values of Storey fc Zeippen ( 2000(1 results in negligi¬ 
ble changes in our Te and abundances calculations. 


Table 5 

Atomic dataset used for recombination lines 


Ion 

Recombination coefficients 

H+ 

Storev & Hummer 119951 

He+ 

Porter et al. 12012. 20131 

He2+ 

Storev & Hummer 119951 

C2+ 

Davev et al. 120001 

0^+ 

Storev 119941 

N2+ 

Fang et al. 12011. 20131 


Table 6 

Percentual contribution of recombination to selected CELs 


Line 

Abell 46 

Abell 63 

Ou5 

[N ii] 5755 

13 

9 

8 

[0 II] 7320,7330 

100 

56 

100 

[0 III] 4363 

14 

43 

34 


22 
























































Table 7 

Ionic abundances and adjs 



Abell 46 

Abell 63 

Ou5 

Ion 


corrected®' 


corrected® 


corrected® 

He+ 

11.15±0.02 

11.15 ±0.02 

11.16±0.02 

11.15 ±0.03 

11.14±0.03 

11.14±0.02 

He2+ 

10.35 ±0.01 

10.35 ±0.01 

10.76 ±0.01 

10.75 ±0.03 

10.78 ±0.01 

10.77 ±0.01 

C2+ (ORLs) 

9.81 ±0.02 

9.81 ±0.02 

9.57 ±0.06 

9.57±0.06 

9.70 ±0.05 

9.70 ±0.05 

N+ 

5.72 ±0.04 

5.78 ±0.04 

6 26 

-0.13 

6.48 

6.32 ±0.05 

6.49 ±0.05 

0+ (CELs) 

6.83 ±0.05 

6.92 ±0.05 

7 1 Q +0.29 
'- 0.22 

7 + 0.22 

‘ -0.18 

6.97 ±0.09 

7.24 ±0.07 

0+ (ORLs) 





8.68 ±0.07 

8.66 ±0.07 

adf{0+) 





51 

26 

0^+ (CELs) 

7.77 ±0.02 

7.85 ±0.02 

017 +u. is 
-0.14 

8.46 

8.05 ±0.06 

8.26 ±0.05 

0^+ (ORLs) 

9.93 ±0.04 

9.93 ±0.04 

Q +0.18 

y.ou _o.20 

9.36 

10.01 ±0.08 

10.01 ±0.08 

adfiO^+) 

145 

120 

15 

8 

91 

56 

Ne2+ 

7.25 ±0.03 

7.34 ±0.03 

7 7 n + 0.22 
-0.18 

8.05 1°;?° 

7.61 ±0.08 

7.87 ±0.05 

S+ 

4.81 ±0.09 

4.86 ±0.08 

4.80 ±0.24 

c; m +0.18 
-0.16 

5.17±0.07 

5.32 ±0.06 

S2+ 

5.90 ±0.05 

5.95 ±0.05 

6.47 ±0.13 

6.63 ±0.13 

6.35 ±0.08 

6.48 ±0.07 

Ar2+ 

5.48 ±0.04 

5.53 ±0.04 

- 0.12 

6 23 +0-^^ 
- 0.11 

5.88 ±0.05 

6.03 ±0.05 

Ar3+ 

5.16 ±0.05 

5.24 ±0.06 



5.72 ±0.08 

5.94 ±0.06 

Ar4+ 





4-48 + 0 ;^^ 

4.66 +o;fo 

C12+ 





4.68 

4.87 + 0 + 

C13+ 





4.57±0.10 

4.72 ±0.10 


Note. —Abundances are indicated in the usual notation as log(X®+/H+)+12. 
After removal of the estimated contribution of recombination to CELs (see text). 
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Table 8 

Total abundances 


Abell 46 Abell 63 Ou5 


Element corrected’^ corrected®' corrected® 


He 

11.21 

± 

0.01 

11.22 

± 

0.01 

11.31 

zb 0.02 

11.30 

zb 0.02 

11.30 

± 0.05 

11.29 

± 

0.05 

N 

6.72 

± 

0.05 

6.77 

± 

0.05 

7.31 

+0.15 

- 0.11 

7.46 

zb 0.11 

7.46 

± 0.05 

7.58 

± 

0.05 

0 

7.85 

± 

0.02 

7.93 

± 

0.02 

8.30 

+ 0.20 

-0.15 

8.59 

+0.16 

- 0.12 

8.18 

± 0.06 

8.40 

± 

0.05 

Ne 

7.29 

± 

0.03 

7.38 

± 

0.03 

7.78 

+ 0.22 

-0.18 

8.14 

+0.19 

-0.15 

7 . 70 : 

TO. 08 
-0.06 

7.96 

± 

0.05 

S 

6.15 

± 

0.05 

6.19 

± 

0.05 

6.75 

zb 0.13 

6.89 

zb 0.12 

6.69 

± 0.08 

6.81 

± 

0.07 

Ar 

5.67 

± 

0.04 

5.72 

± 

0.04 

6.28 

+0.14 

- 0.12 

6.46 

zb 0.12 

6.15 

± 0.05 

6.29 

± 

0.05 


Cl ••• ••• ••• ••• 5.02 5.20 


Note. —All abundances are from CELs except for He. They are indicated as log(X/H)+12. 

® After removal of the estimated contribution of recombination to CELs (see text and Table [Bl). 
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